ABSTRACT
COVID-19 pandemic poses a severe threat to public health. However, so far, there are no effective drugs for COVID-19. Transcriptomic changes and key genes related to Th2 cells in COVID-19 have not been reported. These genes play an important role in host interactions with SARS-COV-2 and may be used as promising target. We analyzed five COVID-19-associated GEO datasets (GSE157103, GSE152641, GSE171110, GSE152418, and GSE179627) using the xCell algorithm and weighted gene co-expression network analysis (WGCNA). Results showed that 5 closely correlated modular genes to COVID-19 and Th2 cell enrichment levels, including purple, blue, pink, tan and turquoise, were intersected with differentially expressed genes (DEGs) and 648 shared genes were obtained. GO and KEGG pathway enrichment analyses revealed that they were enriched in cell proliferation, differentiation, and immune responses after virus infection. The most significantly enriched pathway involved the regulation of viral life cycle. Three key genes, namely CCNB1, BUB1, and UBE2C, may clarify the pathogenesis of COVID-19 associated with Th2 cells. 11 drug candidates were identified that could down-regulate three key genes using the cMAP database and demonstrated strong drugs binding energies aganist the three keygenes using molecular docking methods. BUB1, CCNB1 and UBE2C were identified key genes for COVID-19 and could be promising therapeutic targets.
ABSTRACT
The coronavirus disease 2019 (COVID-19) pandemic causes many morbidity and mortality cases. Despite several developed vaccines and antiviral therapies, some patients experience severe conditions that need intensive care units (ICU); therefore, precision medicine is necessary to predict and treat these patients using novel biomarkers and targeted drugs. In this study, we proposed a multi-level biological network analysis framework to identify key genes via protein-protein interaction (PPI) network analysis as well as survival analysis based on differentially expressed genes (DEGs) in leukocyte transcriptomic profiles, discover novel biomarkers using microRNAs (miRNA) from regulatory network analysis, and provide candidate drugs targeting the key genes using drug-gene interaction network and structural analysis. The results show that upregulated DEGs were mainly enriched in cell division, cell cycle, and innate immune signaling pathways. Downregulated DEGs were primarily concentrated in the cellular response to stress, lysosome, glycosaminoglycan catabolic process, and mature B cell differentiation. Regulatory network analysis revealed that hsa-miR-6792-5p, hsa-let-7b-5p, hsa-miR-34a-5p, hsa-miR-92a-3p, and hsa-miR-146a-5p were predicted biomarkers. CDC25A, GUSB, MYBL2, and SDAD1 were identified as key genes in severe COVID-19. In addition, drug repurposing from drug-gene and drug-protein database searching and molecular docking showed that camptothecin and doxorubicin were candidate drugs interacting with the key genes. In conclusion, multi-level systems biology analysis plays an important role in precision medicine by finding novel biomarkers and targeted drugs based on key gene identification.
ABSTRACT
Severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) infection is a leading cause of pneumonia and death. The aim of this investigation is to identify the key genes in SARS-CoV-2 infection and uncover their potential functions. We downloaded the expression profiling by high throughput sequencing of GSE152075 from the Gene Expression Omnibus database. Normalization of the data from primary SARS-CoV-2 infected samples and negative control samples in the database was conducted using R software. Then, joint analysis of the data was performed. Pathway and Gene ontology (GO) enrichment analyses were performed, and the protein-protein interaction (PPI) network, target gene - miRNA regulatory network, target gene - TF regulatory network of the differentially expressed genes (DEGs) were constructed using Cytoscape software. Identification of diagnostic biomarkers was conducted using receiver operating characteristic (ROC) curve analysis. 994 DEGs (496 up regulated and 498 down regulated genes) were identified. Pathway and GO enrichment analysis showed up and down regulated genes mainly enriched in the NOD-like receptor signaling pathway, Ribosome, response to external biotic stimulus and viral transcription in SARS-CoV-2 infection. Down and up regulated genes were selected to establish the PPI network, modules, target gene - miRNA regulatory network, target gene - TF regulatory network revealed that these genes were involved in adaptive immune system, fluid shear stress and atherosclerosis, influenza A and protein processing in endoplasmic reticulum. In total, ten genes (CBL, ISG15, NEDD4, PML, REL, CTNNB1, ERBB2, JUN, RPS8 and STUB1) were identified as good diagnostic biomarkers. In conclusion, the identified DEGs, hub genes and target genes contribute to the understanding of the molecular mechanisms underlying the advancement of SARS-CoV-2 infection and they may be used as diagnostic and molecular targets for the treatment of patients with SARS-CoV-2 infection in the future.